Linear and nonlinear Stark effect in triangular molecule 
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We analyze changes of the electronic structure of a triangular molecule under the influence of an 
electric field (i.e., the Stark effect). The effects of the fleld are shown to be anisotropic and include 
both a linear and a nonlinear part. For strong electron correlations, we explicitly derive exchange 
couplings in an effective spin Hamiltonian. For some conditions one can find a dark spin state, 
for which one of the spins is decoupled from the others. The model is also applied for studying 
electronic transport through a system of three coherently coupled quantum dots. Since electron 
transfer rates are anisotropic, the current characteristics are anisotropic as well, differing for small 
and large electric field. 
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I. INTRODUCTION 

In this paper we investigate electronic properties of 
a model of a triangular molecule in a presence of an 
electric field. Recently, similar models were considered 
for systems of three coherently coupled quantum dots 
(QDs)?i"— , for magnetic interactions in molecules,--^^ as 
well as for complex phase orderings in strongly corre- 
lated electronic materials with a triangular lattice (e.g., 
multiferroicai^, cobaltatesi^, or organic compounds^-). 
These simple models exhibit plenty of interesting physics. 
For example, in the system of quantum dots, one finds 
complex charge and spin arrangements^, which can be 
classified according a set of topological Hunds rules. 
For some interference conditions a so-called dark state 
can occur, for which one of the QDs is decoupled from 
the reservoirs and an electron can be trappedj^ Con- 
sequently, electronic transport is changed; one can ob- 
serve a rectification effect, negative differential resistance, 
and an enhancement of the shot noise.- Moreover, one 
can expect the Aharonov-Bohm effect, when a magnetic 
flux penetrates the triangle. Moreover, one can expect 
theAharonov-Bohm effect, when a magnetic flux pen- 
etrates the triangle. The effect leads to a crossover 
from the singlet to the triplet ground state, which man- 
ifests itself in spin blockade in transport;^ and interest- 
ing spin dynamics under two-electron-spin-resonance^. 
For the regime of coherent transport one can expect the 
very rich phase diagram with many types of the Kondo 
resonances^. 

Our studies concern the Stark effect in the system 
of strongly correlated electrons and they are addressed 
mainly to coherently coupled quantum dots. For a small 
number of electrons (n = 1 and n = 2 in the triangle) 
the electric field induces a large polarization, and many 
aspects were already considered^"— . In this paper we fo- 
cus on the situation with n — 3 electrons, for which the 
induced polarization is minor, because the Coulomb in- 
teractions dominate and hinder a possible shift of an elec- 
tronic charge. First we will show that the electric field 



leads to splitting of energy levels as well as to breaking 
of the symmetry of the system and changes the symme- 
try of wave functions. We demonstrate that the electric 
field can induce significant changes in a spin arrange- 
ment. It also results in changes of coupling between 
spins and different characteristics of spin-spin correla- 
tion functions with respect to an angle the electric field 
forms with the median of the triangle. In particular we 
will show conditions for the appearance of the dark spin 
state. Second we will show how the Stark effect manifests 
itself in electronic transport. We confine ourselves to the 
case in which at equilibrium the ground state is singlet or 
triplet and for an applied bias voltage excited states with 
three electrons (doublets and quadruplets) participate in 
transport. 



II. INFLUENCE OF ELECTRIC FIELD ON 
SPIN STATES 




FIG. 1: (Color online) The considered model of a triangular 
molecule placed in an electric field E. 



A model of a triangular molecule in the presence of 
an electric field E is shown in Fig.l. The corresponding 



Hamiltonian can be expressed as 



i<j.a 



27r 



where the first term describes electron hopping between 
nearest neighbor sites. For the sake of generality we con- 
sider < < as well as < > 0, for which the model describes 
the system of QDs with electrons or holes, respectively. 
The second and the third terms correspond to insite and 
intersite Coulomb interactions. The last term takes into 
account the influence of the electric field E on the elec- 
tronic polarization P — ^YiUia, where denotes 
the vector pointing to the site e - a charge of an elec- 
tron, 9 - the angle between ri and E, d - the length 
of the arm of the triangle and r — d/VS. One can see 
that the electric field modulates the local site energies: 

= gs cos[9 -I- (i — l)27r/3]. In further considerations we 
take gE = Eer as a coupling parameter. 

For n = 3 electrons in the system the wave functions 
are constructed from the singlet and the triplet states 
by adding an electron (see ilSi). There are quadruplet 
states (^/^QsJ with S = 3/2, = ±3/2, ±1/2, and the 
corresponding wave functions are constructed from the 
triplet states by adding an electron. The ground state 
is, however, the doublet state (^/^DgJ with S = 1/2, 
Sz = ±1/2, which can be formed from the states 

IDsJi = ^cL(cL4 - 44) I vac) ,(2) 



|Dsj2 = ^[2 4-cLcL - cUclAa + 44)] Ivac)(3) 

and the states with double site occupancy c\^c\-ffcl^ | vac) , 

'ia^la'^la l^ac), ctc. The function |D5^)i and |D5^)2 are 
constructed respectively from the singlet and the triplet 
state at the 23-bond by adding an electron to the QD 1. 

Fig. 2 presents an evolution of the spin-spin correlation 
functions when the electric field increases, it manifests a 
transition from the linear to the quadratic Stark effect. 
For small fields the spin-correlation functions oscillate 
with the period 27r. At an intermediate value \gE\ = 4|i| 
a crossover occurs, and for a larger g the functions (Si-Sj) 
show new components due to the quadratic Stark effect. 
One can see in Fig. 2a and 2b that at 6* = (where the 
electric field is perpendicular to the 23-bond and points 
to the 1-st site) the spin correlators (Si • S2) = (S3 • Si) = 
and (S2 • S3) ~ —0.75. This means that the spin at 
QDl is uncoupled from the spins forming the singlet at 
the 23-bond, and the corresponding state is |D5^)i. We 
call such a configuration the dark spin state, in contrast 
to the dark states for n — 1 and n = 2 electrons in the 
triangle molecule, when their properties are connected 
to a specific charge distribution.'*'^ With rotation of the 
electric field, the dark spin state occurs at 9 = 2tt/3 
[9 = 47r/3), when the uncoupled spin is located at the 
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FIG. 2: (Color online) Expectation values of the spin corre- 
lators 81-82, 82 • 83 and 83 • 81 (black, red and blue curves, 
respectively) calculated for the ground state as a function of 
the angle 9 of the electric field with respect to the triangle. 
The plots are obtained for the Hubbard model (1) with a large 
C/o = 30, f/i = 2, i = 1, SB = 0.3 (a), qe = 2.5 (b), qe = ^ 
(c) and = 10 (d). For comparison the dashed curves rep- 
resent the spin correlators for the Heisenberg model (4) with 
the modulated exchange couplings Jij given by Eq.(5) - [for 
smaller Qe < 5 the dashed curves cover the solid ones within 
the plot resolution]. 



QD 2 (3) and the singlet state is at the 12 (31) bond (the 
electric field is then perpendicular to the bond). The 
situation is more complex for a large \gE\ > 4|t|, when the 
quadratic Stark effect begins to dominate and gradually 
changes the period of oscillation of {Si ■ Sj) as well as the 
configuration of the dark spin states. 

In order to understand the crossover from small to 
large fields, we perform the perturbative canonical trans- 
formation of the Hubbard Hamiltonian (1) to an effec- 
tive Heisenberg Hamiltonian, treating the intersite terms 
(both the hopping and the intersite Coulomb interaction 
terms) as small ones.^*" To take into account both the 
linear and non-linear Stark effect the perturbation ex- 
pansion should be carried out to the third order (see the 
Appendix for details). For n — 3 electrons the effective 
Hamiltonian reads 

H = 3Ui+Y, -h, (s, ■ Sj - i") , (4) 

with the exchange coupling 

4i2 4i2(£_g^)2 8i3(2e„_e^_e^.) 



, (5) 



where the indices 1,3,171 denote three different sites. Here, 
we also assumed that the on-site Coulomb interactions 
are stronger than the electric field, i.e. ^ e^. It is 
seen that for the weak field the third order term [the 
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third term in Eq.©] depends linearly on the electric 
field E, whereas the second order term behaves like a 
second power of E and its period of oscillations is twice 
as large as the linear term. This result presents one of 
the main differences between the Stark effect in the sys- 
tem of strongly correlated electrons and that in atomic 
physics^. 

For a very large gE, when the quadratic term domi- 
nates, the dark spin state can occur for E parallel to a 
bond of the triangle. Then the singlet state is formed on 
this bond and the uncoupled spin is at the opposite QD. 
In this limit there is a direct exchange process, which 
is symmetric with respect to exchange of spins between 
QDs, and it does not depend on the orientation of the 
dipole (it is a quadratic dependence on E). Since the 
exchange coupling Jij depends on the difference of site 
energies [see Eq.([5l)], its maximal value is at the elec- 
tric field parallel to the bond. For this case two other 
exchange couplings are equal, and the dark spin state 
appears. 



for which the corresponding wavefunctions are: |D5,)i 
and 1 05^)2 [given by Eq.(l2]) and ([3])], respectively. Fig. 3 
presents the plot of these eigenenergies as a function of 
the electric field. In this case |Ds,)i is the dark spin 
state, but its energy iJoi < £^02 for gs < 4t only. Fig. 3 
also shows the eigenergies Euv and Eu2' for the case 
6 — Tr/3 (when the electric field is perpendicular to the 
13 bond and its direction is opposite to the 2-nd site). 
Now the corresponding eigenstates are 
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However in this case the dark spin state [05^)1/ is the 
excited state {Euv > Eu2')- In the both presented cases 
the states |Ds^)i, (Ds^)!' are constructed from the singlet 
states, whereas |Ds^)2, |Ds^)2' are constructed from the 
triplet states. 




FIG. 3: (Color online) Splitting of states by the electric field in 
the efTective Heisenberg model Q . The solid curves represent 
the states Eui and Eu2 for 9 = 0, whereas the dashed curves 
represent E^v , Eu2' for 9 = tt/3 {t = 1, Uo = 30, Ui = 2). 

The eigenvalues of the Hamiltonian ^ can be readily 
obtained as 



Eq = 3C/i, 
£^Di,2=3C/i-(6J±2A)/4 



for quadruplet 
Here, J = (J12 



and doublet, 
+ J23 + J3i)/3 



(6) 
(7) 

respectively, 
and A = 
~ If the 



V '-^12 + ■^23 + "^31 ~ J12J23 — J12J3I — J23J3I 

electric field is perpendicular to the ij-bond, = e^, 
the exchange couplings J^^ 
spin state. In particular for 9 
eigenvalues are 



Jjm and we have the dark 
0, J12 — J31 and the 



3J 

£"01 = 3C/i — — J23, 
Eu2 = 3J7i — - — J12, 



(8) 
(9) 



III. TRANSPORT IN SEQUENTIAL 
TUNNELING REGIME 

Let us now analyze electronic transport in a sequential 
tunneling regime through our system with the left and 
the right electrode connected to the QD 1 and 2, respec- 
tively. In our calculations we need transfer rates from 
the L (R) electrode to the molecule^^ 



where 



1L{R) ^ |(i^3|cl(2)^|i^2)P/(A£:^2^3 



(12) 



denotes the initial state for n electrons, 
is the corresponding energy dif- 



ference, Jl{r) is a net transfer rate through the potential 
barier between the electrode and the molecule, / denotes 
the Fermi distribution function, the chemical potentials 
are taken as /^l = Ep, /j.r = Ep + eV, Ep is the Fermi 
energy and I^ is a bias voltage. Similarly, one can write 
Tu^^y^ for the reverse tunneling process when the elec- 
tron leaves the molecule. 

Since the quadruplet functions |Qs^) are constructed 
from the triplets |Ts,), therefore, the nonzero transfer 
matrix elements are only: 



l(Q±3/2|cl(2)jT±i)| = V2 |(Q±l/2|cl(2)jTo)| 



V3 |(Q±i 



/2|Ci(2y 



^IT 



±1/ 



N2(l)3l 



(13) 



for a = ±1/2. Here, we use [T^^) (the singlet solution 
|S)) as a linear combination of triplets (singlets) localized 
on the ij bonds, and xjj {xfj) denote the corresponding 
coefficients. The doublet wavefunction is a linear com- 
bination of [Dsji and |Ds^)2 [with the coefficient 
and well as the states with double site occupancy. 

We have checked that for a large Uq the transfer matrix 
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elements for the states with double occupied sites play 
a minor role in electronic transport, and thus they are 
ignored. The corresponding elements are: 



l(D±i/2|cL|T±i)| = V2 |(D±i/2|cL|To)| 

l(D±i/2|cL|T±i)| = V2 |(D±i/2|cL|To)| 
^\{x'^/^/2 + x^/VQ)xI,\, 
|(D±i/2|cL|S)| « Ixf xfal, 

l(D±i/2|cL|S)|«l/2|(xf-V3a:f) 



131 



(14) 

(15) 
(16) 
(17) 

,18 



Next, we solve the corresponding master equatio 

'^^"^ - V p + V r'^- p 

- P-^.i E rt-^.-i + E rt:^.„+J (18) 

to find the occupation probability P^,^ of the eigenstates 
{un) in the steady limit, i.e. for dP^^/dt — 0. The current 
in this limit reads 



p., -r 



(19) 



current 




FIG. 4: (Color online) Maps of the current as a function of the 
bias voltage V and the angle 9 of the electric field with respect 
to the triangle for gs = 0.1 (left column) and qe = 0.8 (right 
column). The Fermi level is taken Ef ~ 5.1, so at equilibrium 
the system contains two electrons in the singlet (triplet) state 
for t = — 1 {t = 1) - top (bottom) plots, respectively. We also 
assumed that Uo = 30, Ui — 2, the temperature T = 0.05, 
7l = 7fl = 70, and the current is in units of 670. 

The results of numerical calculations are presented in 
FigU as a map in the V-6 space. Here we assumed that 



the electric field parameter qe can be controlled indepen- 
dently of the bias voltage, e.g. by means of application 
of additional lateral or back gate electrodes. For small 
fields {gE — 0.1) one can see negative differential resis- 
tance (NDR) (an increase and next a drop of / with V) at 
6* ss and 9 ^ tt. The NDR effect is due to charge accu- 
mulation on the dark state and the interchannel Coulomb 
blockade^. We have also analyzed all current contribu- 
tions 1,^2, 1^3 through various energy levels. As expected 
the current flows via the singlet and the both doublet 
states for t = — 1 (top plots in Fig. 4). Although the 
other states are in the voltage window, they play a minor 
role and the corresponding I^^.i^s exponentially small 
(these processes are thermally activated only) . The situ- 
ation for gE = 0.8 (FigHb) is different. One sees that the 
values of / are higher and the 6 dependence is different. 
At w 7r/3 a new minimum appears. It is an evidence 
of the Stark effect, which is also manifested in activation 
of the triplet and the quadruplet states. The transfer 
of electrons via the quadruplet state is now substantial 
and its contribution /t,q is larger than those from the 
other states. This situation can be explained by ana- 
lyzing the activation energies AE^^^^^ in Fig. 5a and the 
corresponding transfer rates F^jt^j^^. Since the singlet is 
the nondegenerate ground state, it shows the quadratic 
Stark effect and AE'sd is a parabola. In contrast to that 
the activation energies from the triplet state A_Etd and 
AE^Q show, for 9 = Svr/S, a linear and a parabolic de- 
pendence. For this case the triplet levels can be derived 
explicitly as: 

ET = Ui+t-^, 



Et = Ui + -{gs - 2i ± ^9g| + UgEt + 36*2). (20) 

We took the Fermi energy as Ep = 5.1, thus, for a small 
gE, Ai^TD < APtq < Ep and the corresponding trans- 
fer rates F^^j-, and F^^q [see Ea.([T^] are exponentially 
small. For larger fields (gs > 0.2) these energies are 
above Ep and the transfer rates F^l^^ and F^^q are 
activated together with Fg^Q (the current starts to flow 
through all these levels at a threshold voltage). 

The bottom row of Fig. 4 presents the current maps for 
the case t — 1, for which the triplet is the ground state 
at equilibrium. For a bias voltage larger than a thresh- 
old one (for eV ^ 0.7), electrons are transferred via the 
triplet and the doublet states as well as the quadruplet 
state. We can clearly see (at eV w 0.9) the second step 
in the current, when the quadruplet state is activated. In 
this case, the current maps are also different for a small 
and large gE- This results from the dependence of the 
activation energies AEi,^^^^ (and the transfer rates) on 
the electric field. Fig. 5b shows that AEsb (almost) lin- 
early increases with gp and for large gp AEg^j > Ep- 
Therefore, the transfer rate Fgjl^j-, is activated. 

It is worth noticing that the degenerate states (the 
triplet Et for t = —1 as well as the singlet Es for t = 
1) have different dependences on the electric field [see 
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- b) t=\ AE-^^ - 
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FIG. 5: (Color online) Electric field dependence of the energy 
difference AEj^^i^s for the most relevant states participating 
in electronic transport. The plots a) and b) are for the case 
t = —1 and t — 1, which correspond to the current maps in 
Fig. 4 for the top and the bottom panels at ^ = 5tt/3 and 
Ef = 5.1, when at equilibrium the ground state is the singlet 
and the triplet, respectively. Here, we omit the indices 1 and 
2 for the doublet states, because their activation energies are 
very close to each other. 

Ea.pO|. the plots for Atd and Atq in Fig. 5a as well as 
for AsD in Fig. 5b]. One of them is linear vs gE, whereas 
the second one is nonlinear. This is in contrast to the 
Stark effect in atomic physical!, for which all degenerate 
states show a linear field dependence in a wide range. 

Here we analyzed the electronic transport, in which 
two- and three-electron states (with transitions \iy2) ^ 
Ivs)) participated. By using the electron- hole symmetry 
of the model (1) one obtains the same results for tran- 
sitions Ivj) ^ \i'4) (between the states with three and 
four electrons) - provided that one changes the sign of 
the hopping t. 



IV. CONCLUSIONS 

Summarizing, we considered the influence of the elec- 
tric field on strongly correlated electrons in the triangular 
molecule (the linear and the quadratic Stark effect). The 
orientation of E with respect to the molecule is im- 
portant, because the electric field breaks the symmetry 
of the system and changes the symmetry of wave func- 
tions. For some 9 one finds the dark states, responsible 
for negative differential resistance. For n — 3 electrons 
we derived the kinetic exchange coupling J^- between the 
spins, which showed quadratic and linear dependence on 
E. The spin-spin correlation functions exhibit different 
angle 9 characteristics for a small and large E. In partic- 
ular, we studied the dark spin states and their evolution 
with E. The model can be applied to studies of entangle- 
ment of three spin qubitai^ with the electric field in its 
specific role. In particular, it can be applied to a descrip- 
tion of an experiment just recently performed in a system 
of three QDs,^*^ which presented coherent spin manipu- 
lation in a qubit with the logical basis formed from the 



doublet states: ID5J1 and |D5j2- Moreover, we pre- 
dict that the anisotropic Stark effect should be seen in 
electronic transport. 

The studied model is general and can be applied to 
real molecules with the triangular symmetry, to study 
their magnetic and optical features of interest for molec- 
ular spintronics. For example, we predict that spacial 
anisotropy induced by the electric field in the effective 
Heisenberg model will be manifested in molecular mag- 
netism (e.g. in magnetization, magnetic susceptibility, 
or ESR spectra)^ i^°'^^ . Moreover, the model can be the 
paradigm for materials with strongly correlated electrons 
on triangular lattices^"—. In our opinion multiferroics 
are the best candidates to observe the Stark effect, be- 
cause in such materials local ferroelectric orderings can 
modify exchange couplings between spins as well as mag- 
netic orderings. 
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Appendix: Canonical transformation of the Hubbard 
model with modulation of site energy 

Here we apply a canonical perturbation theory for the 
Hamiltonian ([ij for Ui — 0, which we rewrite here as: 

H = W + n+T^i+T+i (A.l) 

where W represent the sum of all the single site terms, 
and T„ is a contribution to the hopping part of the Hamil- 
tonian which increases by n the number of the double 
occupied sites in the system. Using Hubbard operatora^^ 
X"" = |a)(a'| (defined in terms of the exact eigenstates 
I a) for an isolated site i) the contributions to the model 
(lA.ll) can be represented by: 

i.a 

To = J2 UA^fxf^x^xf), 

i<j,a 

T+i = ^ aU.Xf'-'X^-, T_i = Tl (A.2) 

i<j,a 

The site energies: Eia & {0,ei,2ei + Uo} correspond to 
the eigenstates a € {0,t,i, 2}. The following analysis is 
valid for a limit of the small T„, when the differences be- 
tween energies Eia for single and double occupied sites 
are much larger than the hopping parameters. In the 
derivation of the effective Hamiltonian we apply the per- 
turbation theory with respect to the hopping part, which 



6 



can be formulated in terms of the recursive canonical 
transformation.^^ With the use of the Hubbard opera- 
tors the method can be easily generalized to an arbitrary 
form of the on-site zero-order Hamiltonian, including the 
site-dependent terms 

1. The effective second order spin Hamiltonian 

Up to the second order with respect to the hopping 
part the transformed Hamiltonian reads 

H ^e-'^He'^ KHo + ]^[iS,Hx]. (A.3) 

where Ha + Tq and Hi = T+i + r_i. Here, it is 

assumed that in the expansion the linear term with re- 
spect to the off-diagonal part of hopping vanishes, which 
is guaranteed by a condition: 

[iS,Ho\=-Hi. (A.4) 

From this condition we can derive in the explicit form 
the transformation matrix 

iS — ^ ^ iSij , 

_^^o.^2-. ^x^-'^xA . (A.5) 

where A^j = + Uq — tj. After inserting the operator 
5 (jA.SP into (jA.3P we obtain the effective Hamiltonian 
valid up to second order perturbation with respect to the 
hopping part. In a form projected to the subspace Cqq, 
defined as subspace of many-electron states with all sites 
singly occupied, i.e. for n = 3 electrons in the triangle, 
the Hamiltonian reads: 



where 

J^)=2t5(A^i+A^i). (A.7) 

The effective Hamiltonian can be rewritten in a more 
familiar form with a help of the spin operators 

^Icoo -W + Y^jg^ U ■ S, - . (A.8) 



2. The effective third order spin Hamiltonian 

A derivation of the higher order terms in a general case 
is based on the recursive procedure;^ however it is rather 
involved and will be discussed in a separate paper. Here 
we only present the extra 3rd order term projected to the 
Coo subspace 

The exchange parameter reads: 

-hA-^A-i -A-.^A-i - A-.^A-i. - A^^A^M (A.IO) 

' im jm jt mi tj mj mi mj / \ J 

Here, the indices z, j, m denote three different sites. Note, 
that the extra term vanishes for the uniform case ei = 
£2 = £3- 

By expanding jg^ and j[f [Eqs.(A.7) and (A.IO)] in 
a series with respect to e.^ <C Uq we ontain the exchange 
coupling Eq.(5). 
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